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the interchain coupling. For a > a^p, the scaling law is proven to change its functional form, thus 
signaling, for the first time numerically, the onset of coherent transverse two-particle hopping. 
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The physical nature of a system of coupled chains of 
strongly-correlated fermions is currently a very contro- 
versial issue. Such a problem has motivated lots of efforts 
in the recent past, both theoretically and experimentally, 
for a number of fundamental reasons. First, a better 
knowledge of this system will provide further insights to 
understand the dimensional cross-over from one dimen- 
sion (ID) to two dimensions {2D) Secondly, strictly 
ID chains have a very peculiar generic physical behavior 
known as the Luttinger Liquid (LL) behavior and it is es- 
sential to know how stable the LL is with respect to small 
perturbations such as the interchain hopping. Moreover, 
it is not clear yet under which experimental conditions 
the LL behavior can be observed experimentally. 

Some time ago, Anderson suggested |^ that the effect 
of the interchain hopping may be strongly affected by the 
ID character of each chain. It was conjectured that an 
intrachain repulsion of intermediate strength might be 
sufficient to lead to a confinement of the particles within 
each chain. Anderson's confinement scenario has received 
much interest since such a mechanism could explain the 
anomalous transverse transport observed for instance 
in quasi one-dimensional compounds such as the organic 
superconductors 1^,1). 

The LL generic behavior of a ID interacting elec- 
trons chain differs radically from that of a Fermi 
liquid. First, there are no quasiparticle-like excitations 
but rather collective modes with different velocities for 
spin and charge (spin-charge separation). This leads to 
the absence of a step in the momentum distribution at 
the Fermi level but rather to a singularity of the form 
n{k) — n{kp) ~ |fc — kp\°' sign(fci? — k). It is remarkable 



that the exponent a is the only parameter which deter- 
mines completely the low-energy properties of a spinless 
LL. In particular, all the exponents of the static and 
dynamical correlation functions are simply related to a 
(with given sign of the interaction). We shall then con- 
sider a as the key parameter fully determining the im- 
portant properties of the ID metallic system. 

The central issue we shall focus on in the following 
study is the physical role of a small interchain hopping 
t±. Such a question has been addressed by several au- 
thors using different methods and various concepts have 
emerged from these studies such as the notion of rel- 
evance/irrelevance in the Renormalization Group (RG) 
sense or the concept of coherence/incoherence. 

Simple RG calculations suggest that the trans- 
verse hopping is a relevant perturbation for a < 1. In 
that case, the system flows towards a strong-coupling 
fixed point which can not be determined. On the other 
hand, for a > 1, the hopping becomes irrelevant and can 
in principle be neglected. This approach, however, has 
some limitations. First, it is a perturbative method lim- 
ited to first order in t± and there is no guaranty that this 
should work for such a problem. Secondly, even when 
irrelevant, the hopping term always generates new and 
relevant interchain two-particles hopping for all values of 
the LL parameter a ||,^. As a consequence, the system 
always flows to strong coupling and, thus, it seems haz- 
ardous to make predictions about the true ground state 
based only upon the RG arguments. 

Another approach to this problem takes advantage of a 
mapping of the two-chain system onto a two-level system 
coupled to a bath of oscillators jl^ . This study suggests 
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that relevance itself is not a sufficient condition to cause 
coherent motion between the chains. The notion of co- 
herence has been explained in simple terms by Anderson 
and coworkers |]l0| by assuming a system of two separate 
chains prepared at time t — with a different number 
of particles. Then, if the interchain hopping is switched 
on, one can consider the probability P{t) of the system 
returning to its initial state after a time t. Coherence or 
incoherence can then be simply defined as the presence or 
absence of oscillations in P{t). This treatment suggests 
the existence of two different regimes: for a < ag, where 
ao depends on t± and is always smaller than 1/2, co- 
herent motion between the chains takes place while this 
motion becomes incoherent for a > aQ. It is argued 
that, since the interchain hopping is treated as a pertur- 
bation, this result can be applied to an arbitrary number 
of chains. These ideas have been tested extensively by 
numerical methods | [ll| , [T2| showing that the amplitudes 
of the oscillations of P{t) can be drastically affected by 
ergodic properties of the single chain Hamiltonian while 
only the characteristic frequency of the oscillations is a 
reliable measure of interchain coherence. 

In this paper, the role of the interchain hopping is in- 
vestigated by unbiased numerical methods. Exact Di- 
agonalisations (ED) ||l^ of 2 x L (double chains) sys- 
tems of interacting spinless fermions are performed for a 
large set of parameters t± and several system sizes. The 
"ladder" is the simplest geometry which can capture the 
essential mechanisms of the interchain coherence while 
still being tractable numerically. We focus here on sim- 
ple ground-state expectation values related to the basic 
single-particle transverse (i.e. involving charge motion 
between the two chains) Green's functions. In contrast 
to dynamical correlations such as the transverse optical 
conductivity |jl^ such static quantities enable a conve- 
nient finite-size scaling analysis as shown below. Indeed, 
the scaling behavior obtained can be directly compared 
to the ones predicted by various analytical approaches 
hence providing a test of the validity or range of applica- 
bility of these methods. 

First, in Sec. I, we shall describe the model of coupled 
chains with variable-range intra-chain interaction and in 
Sec. II discuss the properties of a single isolated chain. 
In particular, the fundamental ID correlation exponent a 
is calculated as a function of the intra-chain parameters. 
In Sec. Ill, we shall define the difference between the 
momentum distributions 5n{kp) of the two-chain system 
which coincides with the expectation value of the single- 
particle interchain hopping operator of an electron at mo- 
mentum kp and which is the central physical quantity of 
the present analysis. Predictions for Sn(kp) based on var- 
ious analytical approaches will be discussed. In Sec. IV, 
the cross-over to coherent two-particle interchain hopping 
is discussed in terms of the RG flow equations. In Sec. 
V, extensive numerical results are presented for Sn{kp) 
and analyzed using some scaling hypothesis. The scaling 
behaviors based on the numerical results are compared 
to existing analytical treatments. The relevance of more 



complicated two-particle operators is investigated. 

I. THE MODEL 

The model of interacting spinless fermions defined on 
a lattice of two coupled chains of length L can be written 
as follows, 

-^^E(4,is,2+H.c.) (1) 
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where f3 labels the chain (/3 = 1,2), j is a rung index 
(j = 1, . . . , L), Cj^fs is the fermionic operator which de- 
stroys one fermion at site j on the chain /3, and V{r) is 
an intra-chain repulsive interaction between two fermions 
at a distance r (the lattice spacing has been set to one). 
Energies are defined in unit of the intra-chain hopping 
amplitude which has been set to 1. In order to mimic 
a screened Coulomb interaction, we choose a repulsive 
interaction of the form V{i) — 2F/(i -I- 1) for i < iq. 
More specifically, we shall consider here the three cases 
io = 1 , 2 or 3 which correspond to an interaction extend- 
ing up to first, second and third nearest neighbors (NN) 
respectively. For example, in the ig = 3 case, a config- 
uration with two fermions sitting on two lattice sites of 
the same chain at a distance 1, 2 or 3 will contribute to a 
diagonal positive energy of V, |F and V/2 respectively. 
Extending the range of the Coulomb interaction to sec- 
ond and third nearest-neighbor is necessary in order to 
obtain larger values of the exponent a. 

Throughout the paper, we have used closed rings (site 
L is connected to site 1) so that the system is invari- 
ant under discrete translations along the chain direction. 
The "ladder" is then defined on a cylinder. Depending on 
the number of sites, particles, etc... periodic or antiperi- 
odic boundary conditions are used in such a way that 
the corresponding non-interacting system corresponds to 
a closed shell configuration, hence minimizing finite-size 
effects. Ground state properties of these clusters are ob- 
tained by standard ED methods [ p^ . 

II. LL PROPERTIES OF A SINGLE CHAIN 

Before understanding the role of the interchain hop- 
ping we shall first characterize the ID models (i.e. 
t±_ = 0) in terms of a Luttingcr-Liquid description (for 
a coraprehcnsive review concerning this section see e.g. 
Ref. p5| , p^ ). In other words, the charge velocity and 
the correlation exponent a (which are the two important 
physical quantities in the case of spinless fermions) are 
determined as a function of the original parameters of 
the models. 

Crudely speaking, a measures the "force" of the intra- 
chain interaction. However, the range of the interaction 
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also plays major role and a increases sharply with zq as 
seen below. Since the set of ID models as previously 
defined are controlled by two parameters, the magnitude 
and the range of the interaction, and since different mod- 
els can be related to the same value of a, we can then 
investigate in the next Sections whether the anomalous 
dimension a alone controls the interchain transport or 
whether non- universal details of the ID system also mat- 
ter. 

Nevertheless, it is important to notice here that some 
care is needed when working at commensurate densities 
and strong repulsion between fermions. Indeed, when 
the repulsion exceeds a critical value the LL metallic 
phase can undergo a transition to an insulating commen- 
surate Charge Density Wave (CDW) state. In fact, due 
to umklapp scattering, this metal-insulator transition oc- 
curs when the value of a reaches a critical value which 
only depends on the filling factor jist . For low commen- 
surability (i.e. filling factor ^ with large q) one has a 
larger value of the critical a and the metallic LL state is 
then stable in a wider range of the interactions. For this 
reason, we shall consider a density of n = 1/4 where a 
can reach a critical value of about 3. However, we believe 
that the results of this paper are generic and not specific 
to such a filling fraction. 
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FIG. 1. Finite-size scaling of the Drude weight (a) and the 
charge velocity v (b) for the ID spinless Hubbard model with 
y = 4 at n = 1/4 with io = 2 (■) and io = 3 (•). These 
quantities follow a clear behavior. 

Let us here follow the lines of Ref. |Q. For various 
rings of size L, physical quantities such as the Drude 



weight 2ttD {D is the charge stiffness), the charge ve- 
locity V and the compressibility are easily calculated by 
ED methods. Rings with up to 36 sites can be han- 
dled at quarter filling using the Lanczos algorithm. The 
finite-size scaling analysis shown in Fig. m reveals that 
the law expected for a ID LL [0 is very well sat- 
isfied. The extrapolations to the thermodynamic limit 
can then be accurately performed. By using the rela- 
tion 271 D = vK, the exponent a= ^{K+l/K~2 
can be eventually obtained. Results are shown in Fig 
a increases with V but remains small when only NN re- 
pulsion is included. On the other hand, values of a as 
large as 1.5 can be achieved with intermediate values of 
V provided the interaction extends up to distance io = 3. 




FIG. 2. Exponent a ys V for n = 1/4 and for NN interac- 
tion (♦) and longer range io ~ 2 (■), io = 3 (•). 

The ED technique supplemented by finite-size analy- 
sis is itself a very accurate method to investigate model 
Eq. and extract the values of a. Moreover, more con- 
trols on the obtained values of a can be performed. For 
example, the finite-size corrections of the ground-state 
energy per site is predicted by Conformal Invariance ar- 
guments [p^ and (assuming that the central charge is 
equal to 1) is completely determined by v. Similarly, the 
compressibility (which can be directly calculated numeri- 
cally) is uniquely related to v, D and a All these con- 
straints are satisfied numerically in the finite-size scaling 
giving even more confidence in the accuracy of the expo- 
nent a. 

Since the numerical value of a will be crucial in 
the scaling analysis of the next Sections, it is impor- 
tant here to test that correct finite-size scaling behav- 
iors can be obtained for some quantities for a single 
chain. Let us e.g. consider the ground-state correlation 
{(l)o\n{kp)n{—kp)\(f>o) where \(j)o) is the ground state of 
the system and n{k) is the distribution at momentum 
k. This quantity corresponds physically to two processes 
which are depicted schematically in Fig. ||. The first di- 



3 



agram involving an exchange of two particles between 
the two Fermi points at kp and —kp can then be de- 
fined by the connected part (n{kp)n{—kF))^ obtained 
by subtracting the (less interesting) disconnected term, 

{n(kF)n{-kF)) (J = {4)a\n{kF)n{-kF)\(t>Q) 

- {Un{kF)\<i>o){U<-kFMo) ■ (2) 




FIG. 3. Schematic picture of 

the two processes contained in the ground-state correlation 
function < (j)o\n{kF)n{~kF)\4>o >. 

Roughly, one can estimate the large-L behavior of 
this quantity by the following scaling argument: in the 
LL theory, the momentum distribution satisfies n(k) — 
n{kF) ^ \k — /cfl" sign(fci? — fc). However, in a finite 
system of size L the Fermi momentum kp^L) is not pre- 
cisely determined, having an uncertainty of order 1/L 
because of the discreteness of the lattice. Therefore, 
\kF — ~ 1/L and this gives for (n{kF)n{—kF))^ 

a behavior like L^^". We have checked numerically this 
behavior for various models by using the extrapolated 
values of a. For convenience, L(n{kF)n{—kF)') ^ is plot- 
ted in Fig. ^ and shows a very accurate linear behavior 
as a function of L^~^". 




FIG. 4. L < ct)Q\n{kF)n{-kF)\(l)o >c vs L^~'^°' for io = 3, 
V = 1,2, 3, 3.5 from right to left and at density n — 1/4. 



III. SINGLE-PARTICLE TRANSVERSE 
HOPPING 

Transport properties between the chains can be studied 
numerically by considering dynamical correlation func- 
tion such as the transverse Green's function |Q or the 
optical conductivity Although very useful, the nu- 
merical analysis of dynamical correlations is rather in- 
volved and an accurate finite-size scaling is difficult to 
carry out. Here, we shall rather concentrate on ground- 
state equal-time correlations which, while also giving di- 
rect informations on transverse transport, are easier to 
analyse in terms of finite-size scaling. 

A particularly useful physical quantity in the following 
analysis is the momentum distribution n{k,k±) defined 
as usual by 

n{k, ki_) = {MclkAk^ l^t^o) ' (3) 

where the fermion operators c are expressed in the mo- 
mentum representation both for the longitudinal and for 
the transverse momenta. In the case of two coupled 
chains (ladder) the transverse momentum can take the 
two values fcj^ = or fcj^ = tt corresponding to bonding 
or antibonding states. The effect of a small transverse 
hopping tj_ can then be analyzed by considering the dif- 
ference, 

dn{kF) — n{kF,0) — n{kF,T^) (4) 

where kp is the ID Fermi momentum of the t_L = sys- 
tem. The physical meaning of 6n{kF) is clear; it describes 
a single-particle hopping from one chain to the next and 
can also be written as, 

5n{kF) = ('^o|(4^,iC/cp,2 + cL,2CfeF,i)l<^o) : (5) 

where Ck^fi is a destruction operator of a fermion on chain 
(3 with a longitudinal momentum k. 

Since n{k, k±) is simply related to the fermion Green's 
function G{k, kj_, co) by an integration over frequency, we 
expect that n(k,k±) can also give informations on in- 
terchain coherence or incoherence. We shall first briefiy 
review some of the simplest analytical approaches avail- 
able in the literature. As we shall see, different behaviors 
as a function of t± are predicted for n{k, k±) from these 
analytical approaches. Therefore, from a direct compar- 
ison with the analytical behaviors, a numerical analysis 
of n{k, k±) is expected to give useful insights on the rel- 
evant physical mechanisms describing transport of parti- 
cles across the chains. 

One of the simplest perturbative treatment in t±_ which 
has been proposed by Wen |^ and others |^,^|9| con- 
sists in expanding the self-energy in powers of t±. By 
keeping only the lowest order term S(A:, k±,uj) = t±{k±), 
where t±{k±) ~ t±cos{k±), and by using the Dyson 
equation, the Green's function can be written as 
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where Gmik, uj) is the exact Green's function for the iso- 
lated (but fully interacting) single chain. This (RPA-like) 
approximation was shown to become exact for a system 
of an infinite number of chains where each chain is cou- 
pled to all the others |^^. Moreover, the approximation 
is particularly appealing since it gives the two correct 
limits: (i) a = 0, i.e. the formula (^) gives the exact 
free-electron propagator and (ii) t± = 0, i.e. (^) becomes 
the exact ID propagator itself. 

This RPA treatment can be applied to the calculation 
of the momentum distribution difference 5n{kp) in the 
specific case of two coupled chains, i.e. t^{k±^ = 0) = tj. 
and t±{k± = tt) = —t±. The Green's functions for the 
bonding and antibonding states are then given by 



G±(fc,u;) 



1 



[Giz,(fc,c^)]-i±ij 



(7) 



At this step, the form of the \D Green's function is 
needed. For the Tomonaga-Luttinger model with a lin- 
earized dispersion, one can compute this quantity in real 
space ]l6[ . In the case of spinless fermions, the Fourier 
transform can be performed |2l|,^ . Given the fact that 
we are only interested here in dimensional analysis which 
is governed by the anomalous exponent a, we shall take 
the simplest form of the ID Green's function for the right 
movers only: 



G^^{k,uj) cx 



to - 



vk 



(8) 



This function has a branch cut on the real axis for < 
\vk\ {k is defined by k—kp and lo is measured with respect 
to the chemical potential 

Special care is needed for analyzing the analytic prop- 
erties of this Green's function. Introducing a positive 
infinitesimal imaginary part 5 one gets. 



ImG^^ik 



0,LL! + iS) cx 



ReG^^{k^0,uj + i5) cx sgn{uj)\uj\^-" cos{^a) 

(9) 

It follows immediately that the spectral function of the 
two-chain system (proportional to the imaginary part 
of its Green's function) can be written in the form 
A{u;,t±ikj_)) = \u;\"-'^a{t±{k_L)\uj\"-'^sgn{Lj)), where a 
is an a-dependent function. The II? Green's function 
(^ diverges at small frequency when a < 1. In this case, 
a pole in the two-chain Green's function is produced for 
an arbitrarily small t±. A Fermi liquid-like behavior is 
thus recovered with a quasiparticle residue behaving like 

The location of the new poles (measured with respect 
to the chemical potential of the isolated chain) is given 
by the solution of the equation G^^{k,uj = 0) = 0, which 
leads, under the RPA approximation (||), to two real so- 
lutions (one for each sign) for the momentum k corre- 
sponding to two Fermi points kp^ and kp-. This can 



be interpreted as a splitting between the bonding and 
antibonding branches which thus become separated in 
momentum space by 6kp = \kp^ — kp^\. Using the 

previous Green's function, one gets Skp = 1^"^ ""^/v, 
i.e. the Fermi surface warp depends on the strength of 
the electron-electron interaction. This result is actually 
shown to be valid at all orders in t± for the self-energy, 
provided a Fermi surface exists [ p3| . 

Let us now investigate what are the consequences for 
the key parameter Sn{kp). Since the momentum distri- 
bution is given by the integrated spectral function one 
gets, 



5n{kp) ^ j (^A{uj,t^) - A{Lo,-ti^^du 



(10) 



where A is some cut-off proportional to the bandwidth 
or to t (set to 1 for convenience). By introducing a new 

variable of integration x such that w = xt^J_^'^ , we can 
determine the behavior of the integral for small t±^ . When 
a < 1/2, 5n(kp) becomes proportional to times 
a dimensionless integral which is convergent both at small 
and high frequencies so that we can let A — > cx3. However, 
for a > 1/2 a finite cut-off is required to avoid ultraviolet 
divergences and thus, it is found that the dominant term 
in Sn(kp) is linear with t±. 

It is important to stress here that although there exist 
two distinct regimes of scaling of dn{kp) as a function 
of t±_ (namely, for a smaller or larger than 1/2), in both 
regimes there are always real poles in the Green's func- 
tion at two new Fermi momenta away from kp. The 
behavior of 5n{kp) with the interchain hopping t±_ (ob- 
tained for example by numerical methods) is an impor- 
tant quantity giving useful informations on the coupled- 
chain system. Also, it is interesting to note that the be- 
havior 5n{kp) cx predicted by the RPA approach 
when a < 1/2 can also be simply obtained assuming a 
crude picture of two rigid LL momentum distributions 
separated in fc-space by 5k p. Using the well-known re- 
sult for the momentum distribution of a ID LL ||l6|l one 
obtains for small 5k p and for any value of a 



5n{kp) ^ A{5kpY + B5kp 



(11) 



where A and B are a-dependent constant whose expres- 
sion is known. By using the scaling form 5kp ex. t^^^ 
valid in the RPA treatment, and considering that linear 
corrections (not included in this consideration) dominate 
for a > 1/2, one obtains the correct behavior of Sn{kp). 
Of course, this derivation is not completely correct since 
the existence of new Fermi momenta implies that the LL 
form of the momentum distribution is no longer valid 
once t± is finite. 



A different approach has been followed in Ref. |24| by 
calculating directly the linear response to the interchain 
hopping t± of the momentum distribution of an array of 
chains. The main result is the following. 
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n{k, ki_) = nioik) + ti_ cos(fci) (a + B\k - fcFp""^) 

(12) 

for a < 1, and 
n{k,k±) = nioik) +t±cos{kj_)(A + B\k - kpl) (13) 




for a > 1, where nnj^k) is the exact ID momentum 
distribution. If a < 1/2 this formula diverges when k 
approaches kp which indicates the failure of the linear- 
behavior hypothesis at k = kp. On the contrary, for 
a > 1/2, n{k,kj_) and Sn{k) are linear in tj_ also at 
k = kp, in agreement with the RPA results. Although, 
strictly speaking Eq. ( |l^ ) can not be used for fc = A:f 
and a < 1/2, we shall see later that it can nevertheless 
be very useful to interpret our numerical results in the 
t± limit at fixed system length L. 

We finish our brief review by exploring the behavior of 
Sn{kp) within the high-dimensional bosonisation method 
applied to very anisotropic 2D system. It was found |2^] 
that the system of coupled chains is a Fermi liquid with 
a quasiparticle weight Z oc which does not vanish for 
any critical value of a (of course, this is valid only for 
small tj_). The physical picture is very simple consisting 
of two bands separated hy 6kp ^ t±, each band exhibit- 
ing a step-like feature. Therefore, for small t±, the differ- 
ence between the two momentum distributions is directly 
related to the amplitude of the step, 6n{kp) ~ Z. The 
behavior Sn{kp) ~ contrasts with the prediction from 

the RPA 6n{kp) ~ t 



FIG. 5. Non-interacting dispersion relations along the 
chain direction. Open shell configurations for 2 x 16 (left) 
and 2 X 20 (right) clusters. Full (open) symbols correspond 
to occupied (empty) states and fi is the chemical potential. 



dt^/dl = (1 - a)tj_ 
dJ/dl = 2(1 - K)J 



{K -l/K)t]_/2TTVp 



(14) 
(15) 



Using the initial conditions t±{Q) = t± and J(0) = 0, the 
RG flow can be integrated, 



J 



K-l/K 
2'KVp 2a 



2{l-a)l _ 2{1~K)U 



(16) 



a/(l-a) 



needed for further clarifications. 

All previous analytic treatments find, at least for a < 
1, finite quasi-particle residues at some new Fermi points. 
However, one could also wonder whether the effect of the 
transverse hopping could be to generate a splitting be- 
tween the two bands while keeping a LL form. In fact, 
this is indeed the case for some ad-hoc electron-electron 
interaction with equal interchain and intra-chain magni- 
tudes ||2^ , i.e. in which the Fourier transform of the po- 
tential has no component transferring particles from one 
band to the other, or when the chains are connected only 
by density-density interactions and not by hopping 27 1. 



From this expression, the competition between two terms 
can be clearly seen. The first one (associated with 
g2(i-a)i^ is directly related to the one-particle hopping 
while the second term (associated with e^^^"^-*') is re- 
lated to the dimension of the two-particle hopping J . For 
K < a the second term dominates the large-/ (i.e., low- 
energy) behavior and therefore a cross-over is expected 
A numerical study is then when a = a2p = = — 1 ~ 0.41. In section ^ we 



IV. TWO-PARTICLE PROCESSES 

According to the RG analysis applied to this prob- 
lem |^,|| , the single-particle hopping generates under the 
RG flow new processes involving the hopping of two par- 
ticles between neighboring chains: the electron-electron 
(EEPH) and the electron-hole pair hoppings (EHPH). 
The former is relevant for any attractive intra-chain in- 
teraction while the latter becomes relevant for any repul- 
sive interaction. 

As we are interested in the repulsive case, the flows of 
the one-particle hopping t± and of the amplitude of the 
EHPH J are given by the set of coupled equations. 



shall investigate whether this crossover affects the single- 
particle hopping operator. 



V. NUMERICAL RESULTS 

The momentum distribution n{kp,k±) for the two- 
chain model is calculated by diagonalizing exactly by 
means of the Lanczos algorithm a flnite cluster of 2 x L 
sites with L — 8, 12, 16,20 at quarter flUing. Along the 
chains, we use either periodic or antiperiodic bound- 
ary conditions in order to get "open shell" configura- 
tions as defined in Fig. |^. This condition ensures a non- 
degenerate ground state and the possibility of adding or 
removing a particle at the Fermi momentum kp. We 
proceed as follows: first, the absolute ground state of the 
complete Hamiltonian is calculated; then, a new state 
is constructed by applying a destruction operator cor- 
responding to a fermion of momentum {kp,kj_). Even- 
tually, n(kp,k±) is obtained by computing the squared 
norm of the resulting state. The final goal is of course 
to extract an extrapolation to the thermodynamic limit 
from the behavior of Sn{kp) as a function of L. We shall 
see later that such an extrapolation is made possible by 
the existence of a simple scaling function. 

In a finite system of fixed length L we expect to be 
able to write 5n(kp) as a Taylor expansion in powers of 
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t±. Since the change t± —t± leads to the exchange of 
the bonding and antibonding states, this series contains 
only odd powers. For our purpose it is sufBcient here to 
restrict to third order in t± psl, 



Sn{kp) — a{L) t±_ 



h{L)tl . 



(17) 



Here, it is essential to remark that the coefficients a{L) 
and h{L) might depend strongly on the system size. For- 
mally, they can be obtained from the investigation of the 
tj^ limit, e.g. a(L) — '^^g'f^'^ \ t±=Q where the par- 
tial derivative is performed at fixed L. In order to get a 
hint on how dn{kp) should behave in the thermodynamic 
limit, a numerical analysis of the size-dependence of a{L) 
is needed. 

The finite-size dependence of a{L) can, in principle, 
be predicted by applying linear response theory to a fi- 
nite system. In fact, the results of Ref. displayed 
in Eqs. ( ]T^ and (^3|) can be used provided one replaces 
the "cut-off" \k — kp\ with 1/L. Therefore, according to 
Eq. (^2|) and (^3|) , linear response suggests a variation of 
the slope a{L) as for a < 1 and as 1/L for 

a > 1. 
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FIG. 6. Slope a(L) plotted as a function of for 
n — 1/4 and various interactions {io — 3) and for clusters 
of lengths L = 8, 12, 16 and 20. Case a < 1. 

The numerical results for a(L) are shown in Fig. ^ and 
Fig. g for various models and are in perfect agreement 
with the prediction from linear response theory. One can 
notice that different models (i.e with different ranges io) 
with almost the same value of a give very similar results. 
This strongly confirms that only a determines the scal- 
ing law. As seen in Fig. ||, for a < 1/2 the slope a{L) 
diverges with increasing system sizes. One might thus 
expect Snlkp) to vary more rapidly than t± and linear 
response to be no longer valid. The Taylor expansion 
(O) then breaks down in the thermodynamic limit in 
this case. On the other hand, for a > 1/2, a{L) goes to a 
finite limit when L —t oo. It is interesting to notice that 



as one gets close to the cross-over value a = 1 the fits 
in terms of a single power law become less accurate since 
both terms of order 1/L and L^^^" compete with each 
other. 



0.6 



0.4 



0.2 



0.0 



n=1/4 



0.00 



0.05 



0.10 



0.15 



1/L 



FIG. 7. Slope a{L) plotted as a function of 1/L for n = 1/4 
and various interactions and for clusters of lengths L — 8, 12, 
and 16. Case a > 1. 

The fact that a{L) remains finite for increasing L 
(which happens when a > 1/2) is not sufficient to guar- 
anty the validity of the Taylor expansion. To see this we 
have investigated the size dependence of the coefficient 
b{L) of the first non-linear correction. 

The numerical estimations of b{L) show unambigu- 
ously that b{L) dangerously increases with L at least for 
a < 1. Moreover b{L) follows very closely a power law be- 
havior b{L) ^ L^ + const.. The values of the exponents 7 
obtained by a fit of the numerical data are shown in Fig. @ 
and are compared to analytic predictions based on a di- 
agrammatic analysis |2^ . For a smaller than a certain 
value, for which the single-particle hopping is more rele- 
vant than the two-particle one, the previous RPA treat- 
ment (see Sec. Ill, and also the diagrammatic analysis) 
suggests that 7 = 3 — 4q;. In fact, one can show that 
at a given order n in the expansion in t±_ the coefficients 
of the t" term scale like However, when 

two-particle interchain hopping becomes dominant, i.e. 
for a > a2p {a2p = — 1 ~ 0.41), there is a cross-over 
to a different regime. Taking into account the divergence 
at short distances of some diagrams for the self energy 
one obtains 7 = 8 — 2K — 2a ■ Fig. || shows in- 

deed an excellent agreement between these predictions 
and the numerical data. Moreover, it is clear that the 
results do not depend on the details of the model (e.g. 
the range io) but only on the value of a characterizing 
the low-energy behavior. As an additional check, we also 
show in Figs. ^ and |l^ the excellent fits of the numerical 
data with the expected laws in the two regimes. 

This preliminary analysis shows that there exists a 
cross-over to a non-linear regime at large system sizes 
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FIG. 8. Numerical estimations of 7 (obtained from a fit) 
plotted versus a. Models of range io = 3 (•) and io = 2 (□) 
are considered. The analytic predictions (see text) 7 = 3 — 4Qf 
for a < a2p and 7 = 3 — 2a — 2K for a > a2p are shown for 
comparison as full lines. 

(or, equivalently, small temperatures). The cross-over 
takes place when the two terms in 5n{kp) become of the 
same order of magnitude e.g. when ti_ > in the 

regime a < a2p- In the range 1/2 < a < 1, even though 
a{L) has a finite thermodynamic limit, this also happens 
due to the contribution from higher-order diagrams. This 
again signals the instability of the Taylor expansion. 

The failure of the linear response is signaled in Ref. p^ ] 
by the divergence of the linear term at k = kp , whenever 
a < 1/2. This failure occurs, as expected, only at the 
interesting point k = kp. However, even if the coefficient 
a{L) of the linear term has a finite limit, one cannot 
exclude that higher-order terms might become relevant. 
We have indeed shown numerically that this is the case 
for the t\ term. Actually, due to the relevance of t_L, 
higher powers of t±_ carry even more divergent terms in 
the L — > 00 limit. This problem can be also translated 
into the fact that the — > and L — > 00 limits do not 
commute. By considering the linear behavior, we first 
take the tj^ — > limit and we study the size dependence 
of this regime. But, since we are interested in the infinite- 
volume case, we should consider the opposite limit, in 
which L is taken to infinity first, i.e. we shall study 
limfj^^o{limL^oo ^"-(^.f)}- To perform this, we can gain 
some insights from the previous study. Indeed, we have 
obtained the following small-ij^ behavior 

5n{kp) - (ao + aiL^-^^'Y^ - (60 + WL^)t\ , (18) 

where the a.; and hi are L-independent constants. Al- 
though, in principle, ao and 60 could be neglected when 
L ^ 1 and a < 1/2, in practice, for the sizes L we have 
studied and when a gets close to 1/2, it is important to 
consider the oq term in the following analysis. Let us 



FIG. 9. Coefficient b{L) {L = 8, 12, 16 and 20) plotted vs 
^3 -4a ^j^g regime a < a2p. The values of a are shown in 
the plot. 

first focus on the a < a2p regime. Since for a < a2p 
the dominating term at each finite order n of the Taylor 
sum is proportional to one can group these 

elements in terms of a scaling function G'Q[y] of a sin- 
gle variable y = t±L^^°', obtaining for the Taylor sum a 
form GQ[y] L^". Further requiring that the Taylor sum 
has a finite value in the L — > cx3 limit, one can rewrite 
G^y] = G„[y] L"" = tl^^'^'^^ G„[y], or, 

equivalently 

Snikp) - ao <i + t"J^'-"^ F„(L<y(^-")) , (19) 

where the a-dependent function Fa{x) — Ga[x^^°'] (cf. 
Ref. [^) should go to a constant in the a; — > 00 limit, and 
we have restored the linear term that becomes important 
for a close to or larger than 1/2. So far, we have proven 
numerically this scaling form in the regime where the ar- 
gument X of Fa{x) is small, i.e. t± <C L~^^~"\ Indeed, 
using 7 = 3 — 4a, it is easy to rewrite Eq. ( p^ ) into the 
scaling form Eq. (|l9|). If the resummation of the Taylor 
sum as explained above is justified, the scaling form ( p^ ) 
is not restricted to the range a; <C 1 but extends to all 
values of X, in particular to the case x — s- cx) which corre- 
sponds to the thermodynamic limit L — > 00 at fixed t±. 
In this case, as explained above, one expects the func- 
tion Fa to have a finite limit, limj;_^oo Fa{x) = Ca, since 
Snikp) is finite in the infinite-volume limit. Therefore, 

formula (|l9| ) leads naturally to Sn{kp) = Cat"^'^ 
Note that for a < 1/2 the contribution of the aotj. term 
can be neglected for small t± since the exponent a/{l—a) 
is smaller than unity. 

To investigate numerically the validity of the scaling 
relation ( p9| ) for all values of the argument of Fa we pro- 
ceed as follows; from the numerical data 5n{kp) and the 
previous estimations of the constant (which depends 
on a) we construct the quantity 
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in the thermodynamic hmit. 
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FIG. 10. Coefficient b{L) [L = 8, 12, 16 and 20) plotted vs 
j^3-2K-2a ^j^g regime a > a2p. The values of a are shown 
in the plot. 

F'^{L, t^) = {5n{kp) - aotx)A!/^'""^ (20) 
which, a priori is a function of L and <^ independently. 
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FIG. 11. Fa{L,t±) for various values of t_i_ and for lengths 
L = 8, 12, 16 and 20 as a function of a unique variable 

In Fig. |n] is plotted as a function of the combined 
variable LtY^'^ "\ As can be seen on the plot, it is strik- 
ing that, for a < a2p, all the data sets lie on a single 
curve. The scaling hypothesis is then verified to a high 
accuracy. This unique curve then defines the scaling func- 
tion Fa{x) where x = Lt^'"'^ From Fig. [l| it is also 
clear that, when x oo, the function Fa{x) saturates to 
a finite value which, according to the previous discussion, 
implies the asymptotic law 




FIG. 12. Sn{kF) vs t± for various interactions calculated 
on 2 X 8 (o), 2 X 12 (□) and 2 x 16 (0) ladders. The thermo- 
dynamic limit, for a < 1/2, Sn = i'^'^'^ + aotj_ (where 
Ca is estimated from Fig. hll and the linear tj_ term is sub- 
dominant at small t±) is also shown as dashed lines. V = 1, 
2, 3 and 5 correspond to a = 0.10, 0.26, 0.46 and 1.00 re- 
spectively. For V — 5, the dashed line corresponds to a linear 
term only. 

For comparison, we have plotted in Fig. |l^ the raw 
data and the expected L — *■ oo behaviors according to 
( |l9| ) for the quantity Sn{kp) as a function of t±. It is very 
clear from this plot that the finite-size effects are partic- 
ularly strong when t± is small. This result can be qual- 
itatively understood since, as we explained above, the 
thermodynamic limit is obtained when Lt^^^ — > oo. 
A typical length scale Ltyp{t±) is defined from the scal- 
ing behavior and ityp increases rapidly with decreasing 

as 

Our scaling results for a < a2p are in excellent agree- 
ment with the predictions based on the approximate 
(RPA) Green's function Eq. (H), as long as exponents are 
concerned. This agreement is expected to persist at least 
as long as the coherent two-particle interchain hopping 
does not play an important role, i.e. up to the value of 
Q!2p. The approximation first made by Wen that consists 
in neglecting the vertex corrections in the computation 
of the Green's function turns out not to be dramatic as 
proven here numerically, due to the fact that higher-order 
corrections build up in an homogeneous way. 

As stated in Sec. IV, RG calculations [pp| predict a 



(5n(fci?) oc t 



a/(l-a} 



(21) 



cross-over from a one-particle regime to a two-particle 
regime around a — a2p ^ 0.41. Physically, in this 
regime particle-hole hopping dominates with respect to 
single-particle hopping. This cross-over is also signaled 
by the change in behavior of the exponent 7 governing 
the size dependence of b{L) as seen in Fig. |g. In fact, 
in the regime a > a2p, the diagrammatic expansion of 



9 



the self-energy generates nonhomogeneous contributions 
at higher orders in t± so that a re-summation in a sim- 
ple scaling form similar to ( ]l9| ) is quite difficult. By tak- 
ing into account the leading diagrams contributing to the 
self-energy, it has been shown that this crossover changes 
the functional form of the exponent of the behavior of 
Snikp) as a function of t± also in the thermodynamic 
limit. For L 3> 1 and t± <^ t one can argue the scaling 
behavior 

8n{kF) = e/^^-'^^ FKit^l^^^'^'^ L) + aot± (22) 

whose derivation is however not straightforward due 
to the contribution of different inhomogeneous dia- 
grams In this equation, K can be expressed as a 
function of a by inverting the equation a = ^{K+l/K — 
2). In particular, K becomes smaller than a for a > a2p, 
therefore the new exponent a/(l — K) is reduced with 
respect to a/(l — a) and dominates the small-tj^ regime. 
One important consequence of this different scaling be- 
havior is that the anomalous contribution f"^*^^ dom- 
inates with respect to the linear contribution aQt± in a 
larger parameter range, i. e. the behavior of 5n{kp) is 
sublinear up to a = 2/3 (and not only to a = 1/2 as 
obtained within the RPA approximation). This is inter- 
esting since linear response theory, while on the one hand 
predicting its own failure at a < 1/2, due to the diver- 
gence of the coefficient a{L) — uq + aiL^~^" in Eqs. (^ 
[l8| ) , on the other hand would lead to a regular linear be- 
havior for a > 1/2, in contrast with the result of Eq. ([22[). 

Eq. (|2^) has been obtained by cutting the expansion of 
the self energy at a given finite order in t± Due to 
the inhomogeneity of the diagrams, this procedure might 
not produce the correct result, if the Taylor series sums 
up in some unexpected way. It is thus of great importance 
to verify numerically whether there is a deviation at all 
from the scaling behavior Eq. (^9|) for a > a2p and, if this 
is the case, to verify whether the scaling law Eq. (^2|), and 
thus the L ^ oo behavior 5n(kp) = t"^^^ are verified. 
Of course, the deviation of the behavior of the exponent 
7 from the dashed line shown in Fig. ^ already tells us 
that something is changing for a > a2p- However, this 
figure docs not tell us anything about the thermodynamic 
limit. 

The presence of several inhomogeneous contributions 
for a > a2p complicates substantially the numerical anal- 
ysis too. For values of a not too far from a2p (in our 
case for a ~ 0.57) it is difficult to distinguish between 
the two scaling behaviors Eqs. ( p^ and (|2^), due to 
the small difference between the exponents a/{l — K) 
and a/(l — a). Moreover, we shall show that the effects 
of the two-particle contributions start to be dominating 
only at large L, thus forcing us to a careful finite-size 
analysis. In Fig. ([T^), we plot the results of the scal- 
ing for a larger value of a, namely a = 0.7. In curves 
(a) and (b) we proceed in the usual way by plotting 

the quantity = {Sn^kp) — aO tj_)/t'^^^^^^^ as a func- 
tion of — t^^'^ ''■'i, where rj takes the two values 



ry = a in (a) and r] = K in (b), corresponding to the 
two laws Eqs. ( |l9| ) and (|2^), respectively. In both cases, 
the scaling ansatz seems rather poor, thus showing at 
least that something has changed for large a since the 
scaling Eq. ( p^ ) no longer works (Fig. curve (a)). In 
order to improve our accuracy, we further subtract the 
whole linear contribution from Sn^kp) and plot in (c) the 
quantity F'^ j^ — {Sn{kp) — a{L) t±)/t"^^"^~^^ as a func- 
tion of x,j — i^''^ L with rj = K. The subtraction of 
the L-dependent term is harmless in the thermodynamic 
limit, since a{L — > oo) ao for a > 1/2. However, this 
subtraction allows us to eliminate competing terms that 
would make the numerical analysis difficult. This curve 
plotted as (c) in Fig. (|l^) shows that the fit is indeed 
rather good. This shows numerically that for large L the 
scaling Eq. ( ^2|) is the appropriate one. 
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FIG. 13. F^(L,ti_) for 77 = a (a) and ior -q = K (b) and 
Fr],L[L,t±) for r] — K (c) for various values of t± and for 
lengths L = 8, 12, 16 plotted as a function of the variable 
Xrj — tY^^'^^ L (of. text). The interaction V = 4, io = 3 
corresponds to q = 0.7. 

The only flaw of curve (c) in this figure is that it is not 
clear whether Fri{x) goes to a constant in the thermody- 
namic [x — > 00) limit. This should however be expected 
on physical grounds. The reason why this curve does 
not yet saturates is that the system sizes considered are 
still too small to reach the thermodynamic limit in the 
two-particle regime. That is also the reason for which 
it was important to subtract the whole (L-dependent) 
linear contribution to Sn{kp). 

VI. CONCLUSIONS 

In this paper, ground-state correlation functions of 
strongly-correlated coupled chains were investigated by 
numerical exact-diagonalization techniques. First of all, 
the low-energy LL properties of the ID correlated chains 
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were entirely characterized by the Luttinger Liquid cor- 
relation exponent a. The values of a were calculated 
from a finite-size scaling analysis for various strengths 
and ranges of the electron-electron interaction. The cor- 
rect a-dependence of the scaling behaviors of known ID 
correlation functions were recovered. In a second step, 
the expectation value of the single-particle hopping op- 
erator between two coupled chains at k = kp was investi- 
gated by similar ED methods supplemented by finite-size 
scaling analysis. The Taylor expansion of the expectation 
value of the single-particle hopping operator in powers of 
t± was shown to become unstable in the thermodynamic 
limit in agreement with the theoretical prediction that 
the single-particle hopping is relevant. A change of be- 
havior of the size scaling of the coefficient of the t\ term 
for a greater than a critical value a2p is attributed to 
the coherent transverse two-particle hopping becoming 
the dominant perturbation. In addition, in the regime 
a < a2p where transverse two-particle hopping is less rel- 
evant, the finite-size effects can be described in terms of a 
universal scaling function. In the thermodynamic limit, 
it is found that the expectation value of the single-particle 
interchain hopping operator at momentum kp behaves as 



t' 



«/(i- 



in agreement with an RPA-like treatment of the 
interchain coupling. In contrast, in the a > ai-p regime 

a crossover to a i"^*"^ law is observed (dominated by 
a linear contribution when a > 2/3), signaling the dom- 
inance of two-particle hopping processes. 

Whether the coupled-chain system behaves as an or- 
dinary Fermi Liquid is still not clear yet. The energy 
splitting between bonding and antibonding states (which 
should be related to the warping of the Fermi surface) 



calculated numerically in Ref. |14[ varies as t 



.1/(1-") 



as 

suggested by analytic treatments~p3|] . However, for large 
enough a this behavior might occur only above a criti- 
cal value of t_L (see Ref. [Q). Let us also mention that 
transport properties in the direction perpendicular to the 
chains should follow power laws in t j_ . Numerical results 
for the Drude weig ht are indeed compatible with , 
v>2. 
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The single-particle hopping between two chains is investigated by exact-diagonalizations techniques 
supplemented by finite-size scaling analysis. In the case of two coupled strongly-correlated chains 
of spinless fermions, the Taylor expansion of the expectation value of the single-particle interchain 
hopping operator of an electron at momentum kp in powers of the interchain hopping t± is shown 
to become unstable in the thermodynamic limit. In the regime a < a2p {o!2p — 0.41) where 
transverse two-particle hopping is less relevant than single-particle hopping, the finite-size effects 
can be described in terms of a universal scaling function. From this analysis it is found that the 
single-particle transverse hopping behaves as t^^'^ in agreement with a RPA-like treatment of 
the interchain coupling. For a > a2p, the scaling law is proven to change its functional form, thus 
signaling, for the first time numerically, the onset of coherent transverse two-particle hopping. 
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The physical nature of a system of coupled chains of 
strongly-correlated fermions is currently a very contro- 
versial issue. Such a problem has motivated lots of efforts 
in the recent past, both theoretically and experimentally, 
for a number of fundamental reasons. First, a better 
knowledge of this system will provide further insights to 
understand the dimensional cross-over from one dimen- 
sion (ID) to two dimensions {2D) [1]. Secondly, strictly 
ID chains have a very peculiar generic physical behavior 
known as the Luttinger Liquid (LL) behavior and it is es- 
sential to know how stable the LL is with respect to small 
perturbations such as the interchain hopping. Moreover, 
it is not clear yet under which experimental conditions 
the LL behavior can be observed experimentally. 

Some time ago, Anderson suggested [2] that the effect 
of the interchain hopping may be strongly affected by the 
ID character of each chain. It was conjectured that an 
intrachain repulsion of intermediate strength might be 
sufficient to lead to a confinement of the particles within 
each chain. Anderson's confinement scenario has received 
much interest since such a mechanism could explain the 
anomalous transverse transport [3] observed for instance 
in quasi one-dimensional compounds such as the organic 
superconductors [4,5]. 

The LL generic behavior of a ID interacting elec- 
trons chain [6] differs radically from that of a Fermi 
liquid. First, there are no quasiparticle-like excitations 
but rather collective modes with different velocities for 
spin and charge (spin-charge separation). This leads to 
the absence of a step in the momentum distribution at 
the Fermi level but rather to a singularity of the form 
n{k) — n{kp) ~ |fc — kp\" sign(fci? — k). It is remarkable 



that the exponent a is the only parameter which deter- 
mines completely the low-energy properties of a spinless 
LL. In particular, all the exponents of the static and 
dynamical correlation functions are simply related to a 
(with given sign of the interaction) . We shall then con- 
sider a as the key parameter fully determining the im- 
portant properties of the ID metallic system. 

The central issue we shall focus on in the following 
study is the physical role of a small interchain hopping 
Such a question has been addressed by several au- 
thors using different methods and various concepts have 
emerged from these studies such as the notion of rel- 
evance/irrelevance in the Renormalization Group (RG) 
sense or the concept of coherence/incoherence. 

Simple RG calculations [7,8] suggest that the trans- 
verse hopping is a relevant perturbation for a < 1. In 
that case, the system flows towards a strong-coupling 
fixed point which can not be determined. On the other 
hand, for a > 1, the hopping becomes irrelevant and can 
in principle be neglected. This approach, however, has 
some limitations. First, it is a perturbative method lim- 
ited to first order in t±_ and there is no guaranty that this 
should work for such a problem. Secondly, even when 
irrelevant, the hopping term always generates new and 
relevant interchain two-particles hopping for all values of 
the LL parameter a [8,9]. As a consequence, the system 
always flows to strong coupling and, thus, it seems haz- 
ardous to make predictions about the true ground state 
based only upon the RG arguments. 

Another approach to this problem takes advantage of a 
mapping of the two-chain system onto a two-level system 
coupled to a bath of oscillators [10]. This study suggests 
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that relevance itself is not a sufficient condition to cause 
coherent motion between the chains. The notion of co- 
herence has been explained in simple terms by Anderson 
and coworkers [10] by assuming a system of two separate 
chains prepared at time t = with a different number 
of particles. Then, if the interchain hopping is switched 
on, one can consider the probability P{t) of the system 
returning to its initial state after a time t. Coherence or 
incoherence can then be simply defined as the presence or 
absence of oscillations in P{t). This treatment suggests 
the existence of two different regimes: for a < ag, where 
ao depends on t±_ and is always smaller than 1/2, co- 
herent motion between the chains takes place while this 
motion becomes incoherent for a > ag- It is argued 
that, since the interchain hopping is treated as a pertur- 
bation, this result can be applied to an arbitrary number 
of chains. These ideas have been tested extensively by 
numerical methods [11,12] showing that the amplitudes 
of the oscillations of P{t) can be drastically affected by 
ergodic properties of the single chain Hamiltonian while 
only the characteristic frequency of the oscillations is a 
reliable measure of interchain coherence. 

In this paper, the role of the interchain hopping is in- 
vestigated by unbiased numerical methods. Exact Di- 
agonalisations (ED) [13] of 2 x _L (double chains) sys- 
tems of interacting spinless fermions are performed for a 
large set of parameters t j_ and several system sizes. The 
"ladder" is the simplest geometry which can capture the 
essential mechanisms of the interchain coherence while 
still being tractable numerically. We focus here on sim- 
ple ground-state expectation values related to the basic 
single-particle transverse (i.e. involving charge motion 
between the two chains) Green's functions. In contrast 
to dynamical correlations such as the transverse optical 
conductivity [14] such static quantities enable a conve- 
nient finite-size scaling analysis as shown below. Indeed, 
the scaling behavior obtained can be directly compared 
to the ones predicted by various analytical approaches 
hence providing a test of the validity or range of applica- 
bility of these methods. 

First, in Sec. I, we shall describe the model of coupled 
chains with variable-range intra-chain interaction and in 
Sec. II discuss the properties of a single isolated chain. 
In particular, the fundamental ID correlation exponent a 
is calculated as a function of the intra-chain parameters. 
In Sec. Ill, we shall define the difference between the 
momentum distributions Sn{kp) of the two-chain system 
which coincides with the expectation value of the single- 
particle interchain hopping operator of an electron at mo- 
mentum kp and which is the central physical quantity of 
the present analysis. Predictions for Sn{kp) based on var- 
ious analytical approaches will be discussed. In Sec. IV, 
the cross-over to coherent two-particle interchain hopping 
is discussed in terms of the RG fiow equations. In Sec. 
V, extensive numerical results are presented for Sn{kp) 
and analyzed using some scaling hypothesis. The scaling 
behaviors based on the numerical results are compared 
to existing analytical treatments. The relevance of more 



complicated two-particle operators is investigated. 

I. THE MODEL 

The model of interacting spinless fermions defined on 
a lattice of two coupled chains of length L can be written 
as follows, 

- E(4is-,2 + H.c.) (1) 

j 

where [3 labels the chain [[3 = 1,2), j is a rung index 
(j = 1,...,L), Cj^f} is the fermionic operator which de- 
stroys one fermion at site j on the chain [3, and V{r) is 
an intra-chain repulsive interaction between two fermions 
at a distance r (the lattice spacing has been set to one). 
Energies are defined in unit of the intra-chain hopping 
amplitude which has been set to 1. In order to mimic 
a screened Coulomb interaction, we choose a repulsive 
interaction of the form V{i) = 2V / {i -\- 1) for i < ig- 
More specifically, we shall consider here the three cases 
io = 1 , 2 or 3 which correspond to an interaction extend- 
ing up to first, second and third nearest neighbors (NN) 
respectively. For example, in the ig = 3 case, a config- 
uration with two fermions sitting on two lattice sites of 
the same chain at a distance 1, 2 or 3 will contribute to a 
diagonal positive energy of V , and V/2 respectively. 
Extending the range of the Coulomb interaction to sec- 
ond and third nearest-neighbor is necessary in order to 
obtain larger values of the exponent a. 

Throughout the paper, we have used closed rings (site 
L is connected to site 1) so that the system is invari- 
ant under discrete translations along the chain direction. 
The "ladder" is then defined on a cylinder. Depending on 
the number of sites, particles, etc... periodic or antiperi- 
odic boundary conditions are used in such a way that 
the corresponding non-interacting system corresponds to 
a closed shell configuration, hence minimizing finite-size 
effects. Ground state properties of these clusters are ob- 
tained by standard ED methods [13]. 

II. LL PROPERTIES OF A SINGLE CHAIN 

Before understanding the role of the interchain hop- 
ping we shall first characterize the ID models (i.e. 
t±_ = 0) in terms of a Luttinger-Liquid description (for 
a comprehensive review concerning this section see e.g. 
Ref. [15,16]). In other words, the charge velocity and 
the correlation exponent a (which are the two important 
physical quantities in the case of spinless fermions) are 
determined as a function of the original parameters of 
the models. 

Crudely speaking, a measures the "force" of the intra- 
chain interaction. However, the range of the interaction 
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also plays major role and a increases sharply with io as 
seen below. Since the set of \D models as previously 
defined are controlled by two parameters, the magnitude 
and the range of the interaction, and since different mod- 
els can be related to the same value of a, we can then 
investigate in the next Sections whether the anomalous 
dimension a alone controls the interchain transport or 
whether non-universal details of the \D system also mat- 
ter. 

Nevertheless, it is important to notice here that some 
care is needed when working at commensurate densities 
and strong repulsion between fermions. Indeed, when 
the repulsion exceeds a critical value the LL metallic 
phase can undergo a transition to an insulating commen- 
surate Charge Density Wave (CDW) state. In fact, due 
to umklapp scattering, this metal-insulator transition oc- 
curs when the value of a reaches a critical value which 
only depends on the filling factor [15]. For low commen- 
surability (i.e. filling factor ^ with large q) one has a 
larger value of the critical a and the metallic LL state is 
then stable in a wider range of the interactions. For this 
reason, we shall consider a density of n = 1/4 where a 
can reach a critical value of about 3. However, we believe 
that the results of this paper are generic and not specific 
to such a filling fraction. 
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FIG. 1. Finite-size scaling of the Dmde weight (a) and the 
charge velocity v (b) for the ID spinless Hubbard model with 
y = 4 at n = 1/4 with Jo = 2 (■) and lo = 3 (•). These 
quantities follow a clear l/L^ behavior. 

Let us here follow the lines of Ref. [14]. For various 
rings of size L, physical quantities such as the Drude 



weight 27r_D (D is the charge stiffness), the charge ve- 
locity V and the compressibility are easily calculated by 
ED methods. Rings with up to 36 sites can be han- 
dled at quarter filling using the Lanczos algorithm. The 
finite-size scaling analysis shown in Fig. 1 reveals that 
the \jL'^ law expected for a \D LL [17] is very well sat- 
isfied. The extrapolations to the thermodynamic limit 
can then be accurately performed. By using the rela- 
tion [16] 27r_D = vK, the exponent a = \(K -\- 1/K — 2) 
can be eventually obtained. Results are shown in Fig. 2. 
a increases with V but remains small when only NN re- 
pulsion is included. On the other hand, values of a as 
large as 1.5 can be achieved with intermediate values of 
V provided the interaction extends up to distance ig = 3. 
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FIG. 2. Exponent a vs ^ for n = 1/4 and for NN interac- 
tion (♦) and longer range io = 2 (■), io = 3 (•). 

The ED technique supplemented by finite-size analy- 
sis is itself a very accurate method to investigate model 
Eq. (1) and extract the values of a. Moreover, more con- 
trols on the obtained values of a can be performed. For 
example, the finite-size corrections of the ground-state 
energy per site is predicted by Conformal Invariance ar- 
guments [16] and (assuming that the central charge is 
equal to 1) is completely determined by v. Similarly, the 
compressibility (which can be directly calculated numeri- 
cally) is uniquely related to v, D and a [16] All these con- 
straints are satisfied numerically in the finite-size scaling 
giving even more confidence in the accuracy of the expo- 
nent a. 

Since the numerical value of a will be crucial in 
the scaling analysis of the next Sections, it is impor- 
tant here to test that correct finite-size scaling behav- 
iors can be obtained for some quantities for a single 
chain. Let us e.g. consider the ground-state correlation 
(^(j)Q\n(kp)n(—kp)\(j)Q^ where |(/)o) is the ground state of 
the system and n[k) is the distribution at momentum 
k. This quantity corresponds physically to two processes 
which are depicted schematically in Fig. 3. The first di- 
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agram involving an exchange of two particles between 
the two Fermi points at kp and —kp can then be de- 
fined by the connected part (n{kp)n{—kp)')^ obtained 
by subtracting the (less interesting) disconnected term, 

{n{kp)n{-kp))^ = {(j>o\n{kp)n{-kp)\(j>o) 

- {<i>oHkF)\<i>o){Mn{-kF)\h) . (2) 




FIG. 3. Schematic picture of 

the two processes contained in the ground-state correlation 
function < (j)o\n{kF)n{—kF)\(j)o >. 

Roughly, one can estimate the large-L behavior of 
this quantity by the following scaling argument: in the 
LL theory, the momentum distribution satisfies n{k) — 
n{kp) ~ |fc — kp\" sigii{kp — k). However, in a finite 
system of size L the Fermi momentum kp{L) is not pre- 
cisely determined, having an uncertainty of order 1/L 
because of the discreteness of the lattice. Therefore, 
\kp — kp[L)\ ~ 1/L and this gives for (n{kp)n{—kp)')^ 
a behavior like L~^" . We have checked numerically this 
behavior for various models by using the extrapolated 
values of a. For convenience, L(n{kp)n{—kp)') ^ is plot- 
ted in Fig. 4 and shows a very accurate linear behavior 
as a function of L^~'^" . 




FIG. 4. L < (l}o\n{kF)n{-kF)\(l>o >c vs L^'^" for lo = 3, 
V = 1,2, 3, 3.5 from right to left and at density n = 1/4. 



III. SINGLE-PARTICLE TRANSVERSE 
HOPPING 

Transport properties between the chains can be studied 
numerically by considering dynamical correlation func- 
tion such as the transverse Green's function [12] or the 
optical conductivity [14]. Although very useful, the nu- 
merical analysis of dynamical correlations is rather in- 
volved and an accurate finite-size scaling is difficult to 
carry out. Here, we shall rather concentrate on ground- 
state equal-time correlations which, while also giving di- 
rect informations on transverse transport, are easier to 
analyse in terms of finite-size scaling. 

A particularly useful physical quantity in the following 
analysis is the momentum distribution n{k,k±) defined 
as usual by 

n{k,k^_) = {<j>o\clf^^Cf^f^^\<j>o), (3) 

where the fermion operators c are expressed in the mo- 
mentum representation both for the longitudinal and for 
the transverse momenta. In the case of two coupled 
chains (ladder) the transverse momentum can take the 
two values = or = tt corresponding to bonding 
or antibonding states. The effect of a small transverse 
hopping t±_ can then be analyzed by considering the dif- 
ference, 

Sn{kp) = n{kp,0) — n{kp,Tr) (4) 

where kp is the ID Fermi momentum of the t± = sys- 
tem. The physical meaning of Sn{kp) is clear; it describes 
a single-particle hopping from one chain to the next and 
can also be written as, 

Sn(kp) = (<?io|(4^_iCfe^,2 + 4p,2CfeF,i)l'?^o) , (5) 

where is a destruction operator of a fermion on chain 
l3 with a longitudinal momentum k. 

Since n{k, k±) is simply related to the fermion Green's 
function G{k, k±, uj) by an integration over frequency, we 
expect that n{k,k±) can also give informations on in- 
terchain coherence or incoherence. We shall first briefiy 
review some of the simplest analytical approaches avail- 
able in the literature. As we shall see, different behaviors 
as a function of t±_ are predicted for n{k, k±) from these 
analytical approaches. Therefore, from a direct compar- 
ison with the analytical behaviors, a numerical analysis 
of n{k, k±) is expected to give useful insights on the rel- 
evant physical mechanisms describing transport of parti- 
cles across the chains. 

One of the simplest perturbative treatment in t±_ which 
has been proposed by Wen [7] and others [8,18,19] con- 
sists in expanding the self-energy in powers of By 
keeping only the lowest order term T,{k, k±,ui) = t±{k±), 
where t±{k±) = tj_cos(fcj_), and by using the Dyson 
equation, the Green's function can be written as 

G{k,k_i_,ui) = — — - — , , (6) 
[GiD(k,Lj)\ '-+tj_(kj_) 
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where Gioik, uj) is the exact Green's function for the iso- 
lated (but fully interacting) single chain. This (RPA-like) 
approximation was shown to become exact for a system 
of an infinite number of chains where each chain is cou- 
pled to all the others [20]. Moreover, the approximation 
is particularly appealing since it gives the two correct 
limits: (i) a = 0, i.e. the formula (6) gives the exact 
free-electron propagator and (ii) t±_ = 0, i.e. (6) becomes 
the exact ID propagator itself. 

This RPA treatment can be applied to the calculation 
of the momentum distribution difference Sn{kp) in the 
specific case of two coupled chains, i.e. t±{k±_ = 0) = tj_ 
and t±{k±_ = tt) = — The Green's functions for the 
bonding and antibonding states are then given by 



G±(k,ij) 



1 



[GiD{k,L^)]-^±t^ 



(7) 



At this step, the form of the ID Green's function is 
needed. For the Tomonaga-Luttinger model with a lin- 
earized dispersion, one can compute this quantity in real 
space [16]. In the case of spinless fermions, the Fourier 
transform can be performed [21,22]. Given the fact that 
we are only interested here in dimensional analysis which 
is governed by the anomalous exponent a, we shall take 
the simplest form of the ID Green's function for the right 
movers only: 



G^^(fc,w) cx 



(8) 



This function has a branch cut on the real axis for \ui\ < 
\vk\ (k is defined by k — kp and uj is measured with respect 
to the chemical potential //). 

Special care is needed for analyzing the analytic prop- 
erties of this Green's function. Introducing a positive 
infinitesimal imaginary part S one gets. 



ImG^^(k = 0,Lj + iS) cx — l"^! 
ReG~^{k = 0,ui + iS) cx sgn{ui)\ui\'^-" cos{^a) 



w|i-«sin(f a) 



(9) 



It follows immediately that the spectral function of the 
two-chain system (proportional to the imaginary part 
of its Green's function) can be written in the form 
A{ui,t_i_{k_i_)) = \ui\"-'^a{t_i_{k_i_)\ui\"-'^sgn{ui)), where a 
is an a-dependent function. The ID Green's function 
(8) diverges at small frequency when a < 1. In this case, 
a pole in the two-chain Green's function is produced for 
an arbitrarily small A Fermi liquid-like behavior is 
thus recovered with a quasiparticle residue behaving like 

The location of the new poles (measured with respect 
to the chemical potential of the isolated chain) is given 
by the solution of the equation G^^(k,Lj = 0) = 0, which 
leads, under the RPA approximation (6), to two real so- 
lutions (one for each sign) for the momentum k corre- 
sponding to two Fermi points kp^ and kp-. This can 
be interpreted as a splitting between the bonding and 



antibonding branches which thus become separated in 
momentum space hy Skp = \kp^ — kp-\. Using the 



previous Green's function, one gets Skp 



i/(i-«: 



/v, 



i.e. the Fermi surface warp depends on the strength of 
the electron-electron interaction. This result is actually 
shown to be valid at all orders in tj_ for the self-energy, 
provided a Fermi surface exists [23]. 

Let us now investigate what are the consequences for 
the key parameter Sn{kp). Since the momentum distri- 
bution is given by the integrated spectral function one 
gets, 



Sn{kp)= I (^A{Lo,t^)-A{Lo,-t^)^ 



duj 



(10) 



where A is some cut-off proportional to the bandwidth 
or to t (set to 1 for convenience). By introducing a new 
variable of integration x such that uj = xt^J_''^ "\ we can 
determine the behavior of the integral for smalH j_ . When 
a < 1/2, Sn{kp) becomes proportional to t'^^'"^ times 
a dimensionless integral which is convergent both at small 
and high frequencies so that we can let A — cxd. However, 
for a > 1/2 a finite cut-off is required to avoid ultraviolet 
divergences and thus, it is found that the dominant term 
in Sn{kp) is linear with 

It is important to stress here that although there exist 
two distinct regimes of scaling of Sn{kp) as a function 
of tj_ (namely, for a smaller or larger than 1/2), in both 
regimes there are always real poles in the Green's func- 
tion at two new Fermi momenta away from kp. The 
behavior of Sn{kp) with the interchain hopping t±_ (ob- 
tained for example by numerical methods) is an impor- 
tant quantity giving useful informations on the coupled- 
chain system. Also, it is interesting to note that the be- 
havior Sn{kp) cx t^'''"^ predicted by the RPA approach 
when a < 1/2 can also be simply obtained assuming a 
crude picture of two rigid LL momentum distributions 
separated in fc-space by Skp. Using the well-known re- 
sult for the momentum distribution of a ID LL [16] one 
obtains for small Skp and for any value of a 



Sn{kp) =A{Skp)" + BSki 



fll) 



where A and B are a-dependent constant whose expres- 
sion is known. By using the scaling form Skp cx t^J_''^ 
valid in the RPA treatment, and considering that linear 
corrections (not included in this consideration) dominate 
for a > 1/2, one obtains the correct behavior oi Sn[kp). 
Of course, this derivation is not completely correct since 
the existence of new Fermi momenta implies that the LL 
form of the momentum distribution is no longer valid 
once tj_ is finite. 

A different approach has been followed in Ref. [24] by 
calculating directly the linear response to the interchain 
hopping t j_ of the momentum distribution of an array of 
chains. The main result is the following. 
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n(k, kj_) = niD(k) + tj_ cos(fc_L) (a + B\k - fc^P""^) 

(12) 

for a < 1, and 

n{k, kj_) = nioik) + tj_ cos(fc_L) (^A + B\k - kp]^ (13) 

for a > 1, where nio{k) is the exact ID momentum 
distribution. If a < 1/2 this formula diverges when k 
approaches kp which indicates the failure of the linear- 
behavior hypothesis at k = kp. On the contrary, for 
a > 1/2, n[k,ki_) and Sn{k) are linear in t±_ also at 
k = kp, in agreement with the RPA results. Although, 
strictly speaking Eq. (12) can not be used for k = kp 
and a < 1/2, we shall see later that it can nevertheless 
be very useful to interpret our numerical results in the 
tj_ — limit at fixed system length L. 

We finish our brief review by exploring the behavior of 
Sn{kp) within the high-dimensional bosonisation method 
applied to very anisotropic 2D system. It was found [25] 
that the system of coupled chains is a Fermi liquid with 
a quasiparticle weight Z (x t'f_ which does not vanish for 
any critical value of a (of course, this is valid only for 
small tj_). The physical picture is very simple consisting 
of two bands separated hy Skp ^t±, each band exhibit- 
ing a step-like feature. Therefore, for small t±_ , the differ- 
ence between the two momentum distributions is directly 
related to the amplitude of the step, Sn{kp) ~ Z . The 
behavior Sn{kp) ~ tf_ contrasts with the prediction from 

the RPA 5n[kp) ~ t"^''^ "\ A numerical study is then 
needed for further clarifications. 

All previous analytic treatments find, at least for a < 
1, finite quasi-particle residues at some new Fermi points. 
However, one could also wonder whether the effect of the 
transverse hopping could be to generate a splitting be- 
tween the two bands while keeping a LL form. In fact, 
this is indeed the case for some ad-hoc electron-electron 
interaction with equal interchain and intra-chain magni- 
tudes [26], i.e. in which the Fourier transform of the po- 
tential has no component transferring particles from one 
band to the other, or when the chains are connected only 
by density-density interactions and not by hopping [27]. 

IV. TWO-PARTICLE PROCESSES 

According to the RG analysis applied to this prob- 
lem [8,9], the single-particle hopping generates under the 
RG fiow new processes involving the hopping of two par- 
ticles between neighboring chains: the electron-electron 
(EEPH) and the electron-hole pair hoppings (EHPH). 
The former is relevant for any attractive intra-chain in- 
teraction while the latter becomes relevant for any repul- 
sive interaction. 

As we are interested in the repulsive case, the fiows of 
the one-particle hopping t±_ and of the amplitude of the 
EHPH J are given by the set of coupled equations. 




FIG. 5. Non-interacting dispersion relations along the 
chain direction. Open shell configurations for 2 x 16 (left) 
and 2 X 20 (right) clusters. Full (open) symbols correspond 
to occupied (empty) states and jj is the chemical potentiaL 

dtj_/dl=(l-a)tj_ (14) 
dJ/dl = 2(1 - K)J + (K - l/K)i\/2TTVp . (15) 

Using the initial conditions ti_[Q) = tj_ and ■J{Q) = 0, the 
RG fiow can be integrated, 

J^i£-_Vi£^^2(l-«)/_g2(l-K)/)^ (16) 

From this expression, the competition between two terms 
can be clearly seen. The first one (associated with 
g2(i-a);-j jg Jij-gctly related to the one-particle hopping 
while the second term (associated with e^'-"'^""^'") is re- 
lated to the dimension of the two-particle hopping J . For 
K < a the second term dominates the large-/ (i.e., low- 
energy) behavior and therefore a cross-over is expected 
when a = a2p = K2p = — 1 ~ 0.41. In section V we 
shall investigate whether this crossover affects the single- 
particle hopping operator. 

V. NUMERICAL RESULTS 

The momentum distribution n{kp,k±) for the two- 
chain model is calculated by diagonalizing exactly by 
means of the Lanczos algorithm a finite cluster of 2 x _L 
sites with _L = 8, 12, 16, 20 at quarter filling. Along the 
chains, we use either periodic or antiperiodic bound- 
ary conditions in order to get "open shell" configura- 
tions as defined in Fig. 5. This condition ensures a non- 
degenerate ground state and the possibility of adding or 
removing a particle at the Fermi momentum kp. We 
proceed as follows: first, the absolute ground state of the 
complete Hamiltonian is calculated; then, a new state 
is constructed by applying a destruction operator cor- 
responding to a fermion of momentum {kp,k±). Even- 
tually, n[kp,ki_) is obtained by computing the squared 
norm of the resulting state. The final goal is of course 
to extract an extrapolation to the thermodynamic limit 
from the behavior of Sn{kp) as a function of L. We shall 
see later that such an extrapolation is made possible by 
the existence of a simple scaling function. 

In a finite system of fixed length L we expect to be 
able to write Sn{kp) as a Taylor expansion in powers of 
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Since the change t±_ —t±_ leads to the exchange of 
the bonding and antibonding states, this series contains 
only odd powers. For our purpose it is sufficient here to 
restrict to third order in t±_ [28], 

Sn(kF) = a(L)tj_ -b(L)tl . (17) 

Here, it is essential to remark that the coefficients a{L) 
and h{L) might depend strongly on the system size. For- 
mally, they can be obtained from the investigation of the 
tj_ — limit, e.g. a(L) = ^^^^^^\tj_=o where the par- 
tial derivative is performed at fixed L. In order to get a 
hint on how Sn{kp) should behave in the thermodynamic 
limit, a numerical analysis of the size-dependence of a(_L) 
is needed. 

The finite-size dependence of a{L) can, in principle, 
be predicted by applying linear response theory to a fi- 
nite system. In fact, the results of Ref. [24] displayed 
in Eqs. (12) and (13) can be used provided one replaces 
the "cut-off" \k — kp\ with l/L. Therefore, according to 
Eq. (12) and (13), linear response suggests a variation of 
the slope a{L) as {l/L)^"~^ for a < I and as l/L for 
a > 1. 
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FIG. 6. Slope a(L) plotted as a function of for 
n = 1/4 and various interactions (io = 3) and for clusters 
of lengths L = 8, 12, 16 and 20. Case a < 1. 

The numerical results for a{L) are shown in Fig. 6 and 
Fig. 7 for various models and are in perfect agreement 
with the prediction from linear response theory. One can 
notice that different models (i.e with different ranges io) 
with almost the same value of a give very similar results. 
This strongly confirms that only a determines the scal- 
ing law. As seen in Fig. 6, for a < 1/2 the slope a{L) 
diverges with increasing system sizes. One might thus 
expect Sn{kp) to vary more rapidly than t±_ and linear 
response to be no longer valid. The Taylor expansion 
(17) then breaks down in the thermodynamic limit in 
this case. On the other hand, for a > 1/2, a{L) goes to a 
finite limit when L ^ oo. It is interesting to notice that 



as one gets close to the cross-over value a = 1 the fits 
in terms of a single power law become less accurate since 
both terms of order l/L and L^~^" compete with each 
other. 
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FIG. 7. Slope a{L) plotted as a function of l/L for n = 1/4 
and various interactions and for clusters of lengths L = 8, 12, 
and 16. Case a > I. 

The fact that a{L) remains finite for increasing L 
(which happens when a > 1/2) is not sufficient to guar- 
anty the validity of the Taylor expansion. To see this we 
have investigated the size dependence of the coefficient 
h{L) of the first non-linear correction. 

The numerical estimations of h{L) show unambigu- 
ously that h{L) dangerously increases with L at least for 
a < 1. Moreover h{L) follows very closely a power law be- 
havior h{L) ~ L'' -|-const.. The values of the exponents 7 
obtained by a fit of the numerical data are shown in Fig. 8 
and are compared to analytic predictions based on a di- 
agrammatic analysis [23]. For a smaller than a certain 
value, for which the single-particle hopping is more rele- 
vant than the two-particle one, the previous RPA treat- 
ment (see Sec. Ill, and also the diagrammatic analysis) 
suggests that 7 = 3 — 4a. In fact, one can show that 
at a given order n in the expansion in t j_ the coefficients 
of the t1 term scale like [23]. H owever, when 

two-particle interchain hopping becomes dominant, i.e. 
for a > a2p (a2p = V2 — 1 ~ 0.41), there is a cross-over 
to a different regime. Taking into account the divergence 
at short distances of some diagrams for the self energy 
one obtains 7 = 3 — 2K — 2a [9,23]. Fig. 8 shows in- 
deed an excellent agreement between these predictions 
and the numerical data. Moreover, it is clear that the 
results do not depend on the details of the model (e.g. 
the range io) but only on the value of a characterizing 
the low-energy behavior. As an additional check, we also 
show in Figs. 9 and 10 the excellent fits of the numerical 
data with the expected L'' laws in the two regimes. 

This preliminary analysis shows that there exists a 
cross-over to a non-linear regime at large system sizes 
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FIG. 8. Numerical estimations of 7 (obtained from a fit) 
plotted versus a. Models of range io = 3 (•) and io = 2 (□) 
are considered. The analytic predictions (see text) 7 = 3 — 4a 
for a < a2p and 7 = 3 — 2a — 2K for a > a2p are shown for 
comparison as full lines. 

(or, equivalently, small temperatures). The cross-over 
takes place when the two terms in Sn{kp) become of the 
same order of magnitude e.g. when t± > in the 

regime a < a2p. In the range 1/2 < a < 1, even though 
a{L) has a finite thermodynamic limit, this also happens 
due to the contribution from higher-order diagrams. This 
again signals the instability of the Taylor expansion. 

The failure of the linear response is signaled in Ref. [24] 
by the divergence of the linear term at k = kp, whenever 
a < 1/2. This failure occurs, as expected, only at the 
interesting point k = kp. However, even if the coefficient 
a{L) of the linear term has a finite limit, one cannot 
exclude that higher-order terms might become relevant. 
We have indeed shown numerically that this is the case 
for the term. Actually, due to the relevance of 
higher powers of t j_ carry even more divergent terms in 
the L ^ CO limit. This problem can be also translated 
into the fact that the tj_ — and L ^ 00 limits do not 
commute. By considering the linear behavior, we first 
take the tj_ — limit and we study the size dependence 
of this regime. But, since we are interested in the infinite- 
volume case, we should consider the opposite limit, in 
which L is taken to infinity first, i.e. we shall study 
limtj^_>.o{limi_>.oo 5n[kp)}. To perform this, we can gain 
some insights from the previous study. Indeed, we have 
obtained the following small-t j_ behavior 

5n{kp) ~ (ao + aiL^-''")t^ - {bo + bir')tl , (18) 

where the a,- and hi are ^-independent constants. Al- 
though, in principle, ao and &o could be neglected when 
L I and a < 1/2, in practice, for the sizes L we have 
studied and when a gets close to 1/2, it is important to 
consider the ao term in the following analysis. Let us 



FIG. 9. Coefficient h{L) (L = 8, 12, 16 and 20) plotted vs 
^3-4a ^j^g regime a < a2p. The values of a are shown in 
the plot. 

first focus on the a < a2p regime. Since for a < a2p 
the dominating term at each finite order n of the Taylor 
sum is proportional to t'^L"'^^~"^~" , one can group these 
elements in terms of a scaling function Gq,[j/] of a sin- 
gle variable y = t]_L^~" , obtaining for the Taylor sum a 
form Gq,[j/] L~°' . Further requiring that the Taylor sum 
has a finite value in the L ^ 00 limit, one can rewrite 

G„b] = Ga[y] = t"J^'-"^ Ga[y], or, 

equivalently 

Sn{kp) = aot^+ F4Lt'/^'-"^) , (19) 

where the a-dependent function Fa{x) = Ga[x^~"] (cf. 
Ref. [23]) should go to a constant in the x ^ co limit, and 
we have restored the linear term that becomes important 
for a close to or larger than 1/2. So far, we have proven 
numerically this scaling form in the regime where the ar- 
gument X of Fa(x) is small, i.e. t± <^ L~''^~°'\ Indeed, 
using 7 = 3 — 4a, it is easy to rewrite Eq. (18) into the 
scaling form Eq. (19). If the resummation of the Taylor 
sum as explained above is justified, the scaling form (19) 
is not restricted to the range x <^ 1 but extends to all 
values oi X, in particular to the case x ^ 00 which corre- 
sponds to the thermodynamic limit _L — cxd at fixed 
In this explained above, one expects the func- 

tion Fa to have a finite limit, Imix^co Fa{x) = Cq,, since 
Sn{kp) is finite in the infinite-volume limit. Therefore, 

formula (19) leads naturally to 5n[kp) = Cat"^''^ ■ 
Note that for a < 1/2 the contribution of the ao^i term 
can be neglected for smalHj_ since the exponent a /[I — a) 
is smaller than unity. 

To investigate numerically the validity of the scaling 
relation (19) for all values of the argument of Fa we pro- 
ceed as follows; from the numerical data Sn{kp) and the 
previous estimations of the constant ao (which depends 
on a) we construct the quantity 
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in the thermodynamic limit. 
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FIG. 10. Coefficient h{L) (L = 8, 12, 16 and 20) plotted vs 
^3-2A-2a ^j^g regime a > a2p- The values of a are shown 
in the plot. 

F:,{L, t^) = (Snikp) - aot^)/t"J^'-"^ (20) 
which, a prion is a function of L and t±_ independently. 
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FIG. 11. Fa(L, tj_) for various values of tj_ and for lengths 
L = 8, 12, 16 and 20 as a function of a unique variable 
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In Fig. 11 is plotted as a function of the combined 



variable Lt 



1/(1- 



As can be seen on the plot, it is strik- 
ing that, for a < a2p, all the data sets lie on a single 
curve. The scaling hypothesis is then verified to a high 
accuracy. This unique curve then defines the scaling func 
tion Fa[x) where x = Lt 
clear that, when x 
a finite value which, according to the previous discussion 
implies the asymptotic law 



,1/(1- 



j_ . From Fig. 11 it is also 
cxD, the function F„(x) saturates to 



FIG. 12. Sn{kF) vs t± for various interactions calculated 
on 2 X 8 (o), 2 X 12 (□) and 2 x 16 (0) ladders. The thermo- 
dynamic limit, for a < 1/2, Sn = Ca t^^'^ -|- aot± (where 
Ca is estimated from Fig. 11 and the linear t± term is sub- 
dominant at small tj_) is also shown as dashed lines. V = I, 
2, 3 and 5 correspond to a = 0.10, 0.26, 0.46 and 1.00 re- 
spectively. For V = 5, the dashed line corresponds to a linear 
term only. 

For comparison, we have plotted in Fig. 12 the raw 
data and the expected L ^ oo behaviors according to 
(19) for the quantity Sn{kp) as a function of t j_ . It is very 
clear from this plot that the finite-size effects are partic- 
ularly strong when tj_ is small. This result can be qual- 
itatively understood since, as we explained above, the 
thermodynamic limit is obtained when Lt^/'"^ co. 
A typical length scale ityp(^J_) is defined from the scal- 
ing behavior and Ltyp increases rapidly with decreasing 
ti as t 



-1/(1- 



Sn[kp) cx t 



(21) 



Our scaling results for a < a2p are in excellent agree- 
ment with the predictions based on the approximate 
(RPA) Green's function Eq. (6), as long as exponents are 
concerned. This agreement is expected to persist at least 
as long as the coherent two-particle interchain hopping 
does not play an important role, i.e. up to the value of 
a2p. The approximation first made by Wen that consists 
in neglecting the vertex corrections in the computation 
of the Green's function turns out not to be dramatic as 
proven here numerically, due to the fact that higher-order 
corrections build up in an homogeneous way. 

As stated in Sec. IV, RG calculations [8,9] predict a 
cross-over from a one-particle regime to a two-particle 
regime around a = a2p ~ 0.41. Physically, in this 
regime particle-hole hopping dominates with respect to 
single-particle hopping. This cross-over is also signaled 
by the change in behavior of the exponent 7 governing 
the size dependence of h{L) as seen in Fig. 8. In fact, 
in the regime a > a2p, the diagrammatic expansion of 
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the self-energy generates nonhomogeneous contributions 
at higher orders in tj_ so that a re-summation in a sim- 
ple scaling form similar to (19) is quite difficult. By tak- 
ing into account the leading diagrams contributing to the 
self-energy, it has been shown that this crossover changes 
the functional form of the exponent of the behavior of 
Sn{kp) as a function of t±_ also in the thermodynamic 
limit. For L I and t± <^ t one can argue the scaling 
behavior 

Sn(kp) = t"J^'-''^FK(t'/^'-''h) + aot^ (22) 

whose derivation is however not straightforward due 
to the contribution of different inhomogeneous dia- 
grams [23]. In this equation, K can be expressed as a 
function of a by inverting the equation a = ^[K+l/K — 
2). In particular, K becomes smaller than a for a > a2p, 
therefore the new exponent a/(l — A') is reduced with 
respect to a/(l — a) and dominates the small-tj_ regime. 
One important consequence of this different scaling be- 
havior is that the anomalous contribution t^'''"^ dom- 
inates with respect to the linear contribution aot±_ in a 
larger parameter range, i. e. the behavior of Sn{kp) is 
sublinear up to a = 2/3 (and not only to a = 1/2 as 
obtained within the RPA approximation). This is inter- 
esting since linear response theory, while on the one hand 
predicting its own failure at a < 1/2, due to the diver- 
gence of the coefficient a{L) = qq + aiL^~^" in Eqs. (17- 
18), on the other hand would lead to a regular linear be- 
havior for a > 1/2, in contrast with the result ofEq. (22). 

Eq. (22) has been obtained by cutting the expansion of 
the self energy at a given finite order in t±_ [23]. Due to 
the inhomogeneity of the diagrams, this procedure might 
not produce the correct result, if the Taylor series sums 
up in some unexpected way. It is thus of great importance 
to verify numerically whether there is a deviation at all 
from the scaling behavior Eq. (19) for a > a2p and, if this 
is the case, to verify whether the scaling law Eq. (22), and 

thus the L ^ oo behavior 5n[kp) = t"^''^ are verified. 
Of course, the deviation of the behavior of the exponent 
7 from the dashed line shown in Fig. 8 already tells us 
that something is changing for a > a2p. However, this 
figure does not tell us anything about the thermodynamic 
limit. 

The presence of several inhomogeneous contributions 
for a > a2p complicates substantially the numerical anal- 
ysis too. For values of a not too far from a2p (in our 
case for a « 0.57) it is difficult to distinguish between 
the two scaling behaviors Eqs. (19) and (22), due to 
the small difference between the exponents a/(l — A') 
and a/(l — a). Moreover, we shall show that the effects 
of the two-particle contributions start to be dominating 
only at large L, thus forcing us to a careful finite-size 
analysis. In Fig. (13), we plot the results of the scal- 
ing for a larger value of a, namely a = 0.7. In curves 
(a) and (b) we proceed in the usual way by plotting 
the quantity A^ = [5n[kp) — aO t±)/t'^^'"^ '''' as a func- 
tion of = t^/'"^ '^^L, where rj takes the two values 



rj = a in (a) and rj = K in (b), corresponding to the 
two laws Eqs. (19) and (22), respectively. In both cases, 
the scaling ansatz seems rather poor, thus showing at 
least that something has changed for large a since the 
scaling Eq. (19) no longer works (Fig. 13 curve (a)). In 
order to improve our accuracy, we further subtract the 
whole linear contribution from Sn{kp) and plot in (c) the 
quantity A^ ^ = {Sn{kp) — a[L) t\_)lt°'^''^ '''' as a func- 
tion oi x,^ = t^J_''^ '^^ L with T] = K . The subtraction of 
the _L-dependent term is harmless in the thermodynamic 
limit, since a{L co) gq for a > 1/2. However, this 
subtraction allows us to eliminate competing terms that 
would make the numerical analysis difficult. This curve 
plotted as (c) in Fig. (13) shows that the fit is indeed 
rather good. This shows numerically that for large L the 
scaling Eq. (22) is the appropriate one. 
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FIG. 13. Fl,{L,tj_) for r] = a (a) and for r] = K (b) and 
F'^j^(L,tx_) for ri = K (c) for various values of tj_ and for 
lengths L = 8, 12, 16 plotted as a function of the variable 
Xr^ = t^'^ L (cf. text). The interaction V = 4, io = 3 
corresponds to a = 0.7. 

The only fiaw of curve (c) in this figure is that it is not 
clear whether A^ (x) goes to a constant in the thermody- 
namic (x co) limit. This should however be expected 
on physical grounds. The reason why this curve does 
not yet saturates is that the system sizes considered are 
still too small to reach the thermodynamic limit in the 
two-particle regime. That is also the reason for which 
it was important to subtract the whole (_L-dependent) 
linear contribution to Sn{kp). 

VI. CONCLUSIONS 

In this paper, ground-state correlation functions of 
strongly-correlated coupled chains were investigated by 
numerical exact-diagonalization techniques. First of all, 
the low-energy LL properties of the ID correlated chains 



10 



were entirely characterized by the Luttinger Liquid cor- 
relation exponent a. The values of a were calculated 
from a finite-size scaling analysis for various strengths 
and ranges of the electron-electron interaction. The cor- 
rect a-dependence of the scaling behaviors of known \D 
correlation functions were recovered. In a second step, 
the expectation value of the single-particle hopping op- 
erator between two coupled chains at k = kp was investi- 
gated by similar ED methods supplemented by finite-size 
scaling analysis. The Taylor expansion of the expectation 
value of the single-particle hopping operator in powers of 
t±_ was shown to become unstable in the thermodynamic 
limit in agreement with the theoretical prediction that 
the single-particle hopping is relevant. A change of be- 
havior of the size scaling of the coefficient of the term 
for a greater than a critical value a2p is attributed to 
the coherent transverse two-particle hopping becoming 
the dominant perturbation. In addition, in the regime 
a < a2p where transverse two-particle hopping is less rel- 
evant, the finite-size effects can be described in terms of a 
universal scaling function. In the thermodynamic limit, 
it is found that the expectation value of the single-particle 
interchain hopping operator at momentum kp behaves as 
«) agreement with an RPA-like treatment of the 
interchain coupling. In contrast, in the a > a2p regime 

a crossover to a t'^^'"^ law is observed (dominated by 
a linear contribution when a > 2/3), signaling the dom- 
inance of two-particle hopping processes. 

Whether the coupled-chain system behaves as an or- 
dinary Fermi Liquid is still not clear yet. The energy 
splitting between bonding and antibonding states (which 
should be related to the warping of the Fermi surface) 

calculated numerically in Ref. [14] varies as t^/'"^ as 
suggested by analytic treatments [23]. However, for large 
enough a this behavior might occur only above a criti- 
cal value of t j_ (see Ref. [14]). Let us also mention that 
transport properties in the direction perpendicular to the 
chains should follow power laws in Numerical results 
for the Drude weight [14] are indeed compatible with , 
> 2. 
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